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ABSTRACT 


Experimental  data  from  Davidson  Laboratory  on  induced  mass  and  drag 
of  the  basic  finner  missile  have  been  reevaluated  at  the  Naval  Weapons 
Laboratory.  All  but  four  of  the  runs  must  be  discarded  as  determina¬ 
tions  of  induced  mass  because  of  uncertainty  in  velocity.  It  is  con¬ 
cluded  that  an  insight  into  the  mechanism  of  induced  mass  can  only  be 
achieved  thrpugh  mathematical  analysis. 

FOREWORD 

This  report  has  been  prepared  in  compliance  with  BtiWEPS  Directive 
RRRE~07  004/210-1/R009-03-003  dated  10  October  1960.  -Date  of  com¬ 
pletion  was  6  April  1961. 

DfraopocnoN 

An  extensive  effort  has  been  made  in  various  laboratories . to  deter¬ 
mine  the  hydrodynamic  charact<  fistics  of  missiles.  .Quantitative  data 
on  forces  and  moments  are  available  for  steady  flight.  The  interpre¬ 
tation  and  utilization  of  data  for  unsteady  flight  require  a  know¬ 
ledge  of  the  induced  mass  and  moments  of  inertia  of  the  entrained 
fluid.  Theoretical  studies  of  this  induced  mass  have  been  limited  to 
ideal  fluids.  No  reliable  information  is  available  yet  about  how  the 
induced  mass  varies  with  Reynolds  number  or  acceleration  modulus.. 

The  California  Institute  of  Technology  has  investigated  the  motion 
of  the  basic  finner  missile  in  vertical  and  horizontal  flight. 
Quantitative  interpretation  of  the  horizontal  runs  probably  will  be 
possible  after  we  know  the  dependence  of  induced  mass  on  Reynolds 
number.  Quantitative  correlation  of  the  vertical  runs  has  not  been 
possible  without  an  assumption  that  the  actual  masses  were  incor¬ 
rectly  recorded.  No  reruns  have  been  possible  because  the  apparatus 
has  been  dismantled. 

The  Naval  Weapons  Laboratory  proposed  in  reference  1  that  programmed 
acceleration  trials  be  run  in  a  towing  tank.  The  Davidson  Laboratory 
has  made  such  trials  and  has  established  new  values  for  the  steady 
state  drag  of  the  basic  finner  missile.  The  drag  curve  has  a 
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mysterious  hump,  but  no  documentation  of  flow  regimes  is  available.  . 
to  explain  the  cause  of  the  hump.  The  instrumentation  was  barely 
adequate  to  measure  the  induced  mass  ( k ^  **  0. 15) .  The  following 
sources  of  error  are  considered  noceworthy. 

a.  The  interior  of  the  model  was  flooded  with  fluid.  Although 
a  portion  of  the  interior  was  plugged,  the  pressure  gradients  during 
acceleration  were  those  of  a  hollow  shell  filled  with  fluid.  The 
equivalent  mass  of  the  fluid  content  was  a  substantial  part  of  the 
internal  mass,  and  this  caused  a  loss  of  sensitivity. 

b.  The  interior  of  the  model  was  exposed  to  the  fluid  pressure 
at  the  rim  of  the  flat  base.  The  steady  drag  was  correctly  deter¬ 
mined  insofar  as  the  base  pressure  is  uniform.  The  effect  of 
induced,  mass  was  probably  in  error  because  the  pressure  distribution 
from  induced  mass  is  not  uniform  over  the  base. 

c.  The  velocity  of  the  model  could  not  be  controlled  to  follow 
closely  an  ideal  stepwise  variation.  The  velocity  records  were 
reccrded  at  too  small  a  scale  to  be  read  with  precision. 

Although  these  errors  largely  obscure  the  induced  mass,  a  reanalysis 
has  been  made  at  the  Naval  Weapons  Laboratory  in  an  attempt  to 
retrieve  some  useful  information. 

TRAJECTORY  INTEGRATIONS 

The  test  conditions  of  the  runs  were  forwarded  by  Davidson  Laboratory 
in  reference  '2  and  the  original  records  were  forwarded  in  refer-  *. 
ence  ’3.  A  catalog  of  the  test  conditions  is  given  herewith  in 
Table  I  and  samples  of  the  records  are  reproduced  in  Figures  1  and  ’2. 

With  the  exception'  of  the  first  four  runs  in  Table  I  the  velocities 
were  recorded  at  much  too  small  a  scale.  The  smallest  scale  divi¬ 
sion  on  the  speed  records  has  been  converted  to  speed  increments 
through  the  use  of  a  calibration  chart  which  was  supplied  with 
reference  -3.  The  speed  .'ncrement  per  scale  division  is  listed  in 
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Table  I.  Even  with  the  assumption  that  the  speed  records  can  be 
estimated  correctly  to  one  fifth  of  a  scale  division  the  uncertainty 
in  speed  is  enough  to  mask  the  induced  mass*.  Thus  there  would  be  a 
lOOJc  error  in  induced  mass  if  there  were  an  error  of  0.36  (ft)-/(sec) 
over  a  velocity  range  of  5  (ft) /(sec). 

The  uncertainty  in  the  records  from  wiggles  is  clearly  apparent  from 
an  inspection  of  Figures  1  and  2. 

The  original  records  from  the  Davidson  Laboratory  have  been 
reanalyzed  at  the  Naval  Weapons  Laboratory  on  the  basis  of  a 
momentum-displacement  correlation  instead  of  a  force-velocity  corre¬ 
lation.  The  objective  of  the  change  in  correlation  was  a  reduction 
in  the  uncertainty  from  wiggles  in  the  oscillographic  records.  The 
drag  records  were  reevaluated  to  obtain  a  set  of  average  drag  forces 
whose  summation  would  reproduce  the  area  under  the  drag  traces. 

The  results  of  this  reevaluation  are  documented  in  Appendix  C. 

Simplified  trajectory  integrations  were  performed  on  ftORC  with  a 
uniform  time  interval  of  0.5  (sec).  The  trajectory  integrations 
utilized  a  mass  coefficient  m.  and  a  drag  coefficient  k .  The  mass 
coefficient  m  (slugs)  for  the  test  models  is  'related  to  the  added 
mass  coefficient  k1  by  the  equation** 

n  =  (1.01)  +  (0.50)  kt  (l) 

and  the  drag  coefficient  k  (slugs)/ (ft)  is  related  to  the  drag 
coefficient  Cp  by  the  equation** 

k  =  (0.0888)  CD  (-2) 

where  Cp  applies  to  the  base  area  alone. 

*ln  reference  1  the  specification  of  accuracy  was  0.5$  of  full 
scale. 

**These  equations  are  derived  from  a  base  area  of  0.0916  (ft)2,  a 
total  volume  of  0,'2599  (ft)8,  and  a  total  weight  of  -32.48  (lb), 
as  quoted  in  reference  4. 


The  accumulated  error  in  impulse  en  (lb)  (sec)  after  n  steps  of 
integration  is  given  by  the  equation 


where  7t-  (lb)  is  the  average  force  in  the  tth  interval  and  vi  (ft) /(sec) 
is  the  velocity  at  the  end  of  the  ita  interval.  The  coefficients 
were  interpolated  from  a  table  of  values  at  half  intervals  in  velo¬ 
city.  A  table  of  values  is  given  in  Appendix  C.  Thus  the  values 
n-  were  computed  from  the  equation 

mi  =  ny  +  (>2  Vi  -  j)  (my+1  -  7Hj)  (4) 

and  the  values  k-  were  computed  from  the  equation 

ki  =  kj  +  (2  vi  -  j)  Uy+1  -  kj)  (5) 

where  j  is  the  serial  number  in  the  table  for  that  entry  Vj  which 
is  next  smaller  than  v The  results  of  computation  were  plotted  on 
the  CRT  printer.  A  set  of  results  is  given  in  Appendix  D  where  £„ 
is  plotted  against  vn,  and  each  point  is  labeled  with  the  value  of 
n.  An  error  in  mj  is  reflected  in  the  plots  by  a  nonzero  slope 
during  increment  of  velocity  and  an  error  in  kj  is  reflected  by  a 
nonzero  drift  rate  during  stationary  velocity.  The  entries  in  the 
table  of  coefficients  are  so  adjusted  by  trial  as  to  minimize  the 
random  deviation  of  the  plots  from  the  velocity  axis. 

Various  adjustments  of  mass  and  drag  were  tested.  The  induced  mass 
could  be  varied  by  '25$  .from  the  value  which  is  reported  in  refer¬ 
ence  4  without  appreciable  improvement  in  the  error  curves.  The 
drag  from  the  constant  speed  runs  was  not  the  optimum  and  an  improve¬ 
ment  in  the  error  curves  could  be  achieved  through  an  adjustment  of 
drag.  Values  of  kj  from  the  constant  speed  runs  are  listed  in  the 
second  column  of  Table  II  and  the  values  of  kj  after  adjustment  are 
listed  in  the  third  column  of  the  table*.  The  adjusted  values  are 
basic  to  the  error  curves  in  Appendix  D. 

*Although  constant  values  are  listed  in  the  table  for  low  velocity 
and  for  high  velocity  these  were  never  used  in  the  actual  integra¬ 
tions. 


Even  if  all  runs  are  rejected  except  the  first  four  because  of  error 
in  velocity,  there  is  still  a  discrepancy  in  induced  mass  between 
runs  2  and  20. 

FLOW  ANALYSIS 

Experimental  determinations  of  induced  mass  heretofore  .hsive  not  been 
quantitative.  It  seems  obviou*  that  mathematical  analysis  is  ner.es-. 
sary  for  an  insight  into  the  characteristics  of  induced  mass,. 

If  a  missile  were  accelerated  suddenly  from  one  constant  velocity  to 
another  then  a  potential  flow  would  be  superimposed  upon  the  preexisting 
flow.  Since  the  potential  flow  would  not  satisfy  the.  boundary  condi¬ 
tion  of  constant  velocity  at  the  surface  of  the  missile,  .the  potential 
flow  would  be  modified  gradually  through  a  diffusion  of  vorticity. 
Meanwhile  the  drag  would  decay  from  a  large  initial  value  to  a  steady 
final  value.  The  acceleration  thus  would  initiate  a  greater  total 
impulse  than  that  required  to  create  the  potential  flow. 

If  a  rapid  acceleration  cycle  were  applied  to  the  missile  the  diffu¬ 
sion  of  vorticity  would  not  have  time  to  develop  and  the  induced 
mass  would  be  just  that  of  the  classical  potential  flow.  A  com¬ 
puting  program  for  potential  flow  over  missiles  has  been  developed 
by  the  Douglas  Aircraft  Company.  Details  are  given  in  reference  5 
and  subsequent  reports.  The  computing  program  can  be  applied  to  the 
basic  -firmer  missile. 

If  a  slow  acceleration  cycle  were  applied  to  the  missile  the  diffu¬ 
sion  of  vorticity  would  have  time  to  develop  and  the  total  impulse 
would  be  the  integrated  result  of  differential  increments  of  flow 
configuration. 

A  theoretical  study  of  induced  mass  is  underway  at  the  Naval  Weapons 
Laboratory.  An  initial  model  consisted  of  a  pair  of  line  vortices 
behind  a'  cylinder.  The  effect  of  the  vortices  was  found  to  be  a 
decrease  of  induced  mass.  An  acceleration  moves  the  vortices  closer 
to  the  cylinder,  and  diminishes  the  pocket  of  entrained  fluid. 


5 


Although  circulatory  motion  can  be  observed  in  the  trailing  wakes  • 
of  cylinders,  the  concept  that  vorticity  is  concentrated  in  the 
pocket  is  illusory.  Valid  solutions  of  the  Navier-Stokes  equation 
show  that  the  vorticity  trails  off  from  the  cylinder  in  a  vortex 
sheet  from  each  separation  point.  There  is  relatively  little  vor¬ 
ticity  in  the  pocket  of  entrained  fluid,  while  there  is  even  less 
potential  gradient. 

A  new  computing  program  for  solving  the  Navier-Stokes  equation  is 
now  in-  preparation.  The  new  program  will  give  time  dependent  solu¬ 
tions  for  flow  past  a  cylinder.  Line  vortices  are  placed  at  the 
intersections  in  a  grid.  The  rate  of  change  of  vortex  strength  at 
each  grid  point  is  determined  by  finite  difference  approximations  of 
the  diffusion  and  convection  of  vorticity.  Stream  function  is  deter¬ 
mined  by  the  summation  of  contributions  from  each  line  vortex. 

Storage  requirements  in  the  calculator  are  minimized  by  the  use  of  a 
polar  grid.  This  program  will  provide  the  first  determination  of  a 
variation  of  induced  mass  with  Reynolds  number. 

More  information  about  flow  regimes  is  needed.  Possibly  small 
models  of  the  basic  finner  missile  could  be  moved  through  a  bentonite 
suspension  in  a  tank  with  polaroid  windows.  Photographs  of  the" 
double  refraction  would  show  the  onset  of  turbulence  at  various 
points  on  the  jnissile. 

RECOMMENDATIONS 

1.  It  is  recommended  that  further  tests  on  induced  mass  be  spon¬ 
sored  at  Davidson  Laboratory,  but  only  if  all  of  the  following 
specifications  are  met; 

a.  The  power  drive  and  recording  system  be  modified  to  give  a 
better  control  and  a  more  precise  determination  of  the  velocity. 

b.  The  model  be  mounted  on  side  struts  instead  of  the  basT 
sting  (as  recommended  by  Davidson  Laboratory) . 

c.  The  interior  of  the  model  be  sealed  off  from  fluid  contact 
(with  dynamometer  in  struts). 


2_  it  is  reccrmendea  that  a  project  be  established  at  a  hydraulic 
laboratory  for  the  photography  of  flow  regimes. 

3.  It  is  reccrosnded  that  a  project  be  established  at  Douglas 
Aircraft  Company  for  the  computation  of  potential  flow  over  the 
basic  finner  missile. 

4.  It  is  reccnmended  that  the  programming  and  calculation  of  the 
vortex  strength  behind  a  cylinder  be  continued  at  the  Naval  Weapons 
Laboratory  to  the  point  of  determining,  the  induced  mass  of  the 
entrained  fluid. 

CD?  "LgSZCEf 

It  is  concluded  that  an  insight  into  the  mechanism  of  induced  mass 
will  not  be  gained  without  a  mathematical  analysis  of  flow  regimes. 
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APPENDIX  A 


TABLES 


TABLE  I 


CATALOG  OF  RECORDS  FROM  ORIGINAL  BATA 


Run 

Number 

Sting 

Diameter 

(in) 

Speed 

Range 
(ft)-/  (sec) 

Scale 
Division 
(ft)-/  (sec) 

Run 

Number* 

•2 

1.75 

0  - 

11.28 

0.12 

1 

6 

1.75 

0  - 

10.45 

0.12 

•2 

9 

1.75 

0  - 

17.22 

0.12 

-3 

20 

1.75 

6. 12  - 

11. 14 

0.12 

5 

21 

1.75 

11. 76  - 

16.50 

1.8 

6 

22 

1.75 

9.85  - 

5.8S 

1.8  ' 

7 

25 

1.75 

14.06  - 

8.40 

1.8 

9 

2? 

1.75 

18.18  - 

14.00 

1.8 

10- 

•28 

1.75 

9.70  - 

6.50 

1.8 

'•  8 

•33 

■2.50 

O  - 

11.96 

1.8 

13. 

-33 

2.50 

11.96  - 

5.72 

1.8 

13 

39 

2.50 

1-3.63  - 

17.50 

1.8  ■ 

12 

47 

2.50 

17.69  - 

14'.  00 

1.8 

14  • 

54- 

1.75 

0  - 

16.22 

1.8 

-4‘ 

*Runs  as  renumbered  in  the  final  report,  reference  4. 
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TA3LS  II 


DRAG  COSFFICiaTIS  FROM  CCNSSAKT  SPEED  HOIiS  AND 
AFTER  ADJUSUSSrr  FOR  OPTI1SOJI  ERROR  CORTES 


l 

Vi 

k3 

Constant  Speed 

k3 

Optimum  Ei 

- 

0 

0 

0.0372 

0.0372 

•2 

1 

0.0372 

0.0372 

4 

2 

0.0372 

0.0372 

6 

-3 

0.0372 

0.0372 

8 

4* 

0.0363 

0.0363 

10 

5 

0.0354 

0.0354 

12 

6 

0.0346 

0.0346 

14 

7 

0.0341 

0.0341 

16 

8 

0.0381 

0.0381 

18 

9 

0.0392 

0.0392 

•20 

10 

0.0375 

0.0354 

•22 

11 

0.0354 

0.0353 

24 

12. 

0.0352 

0.0352 

26 

13 

0.0347 

0.0347 

28 

14 

0.0344- 

0.0344- 

-30 

15 

0.0341 

0.0344 

-32 

16 

0.0340 

0.0344- 

-34,' 

17 

0.0339 

0.0344- 

36 

18 

0.0339 

0.0344- 

•38 

19 

0. 0338 

0.0344- 

40 

■20 

0.0338 

0.0344 

* 
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Figura  1  -  Suqpla  Spaad  Racorda 
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APPENDIX  D 


Plots  of  en  (ordinates)  versus  vn  (abscissae)  for  each  run 
with  each  point  marked  by  the  value  of* n.  (Sting  diameter  =  1.75 
in.  The  data  for  the  larger  sting  diameter  ran  off  scale. ) 
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